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Abstract 

We apply a simple analytical criterion for locating critical temperatures to Potts 
models on square and triangular lattices. In the self-dual case, i. e. on the square 
lattice we reproduce known exact values of the critical temperature and derive the 
internal energy of the model at the critical point. For the Potts model on the 
triangular lattice we obtain very good numerical estimate of the critical temperature 
and also of the internal energy at the critical point. 
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1 Introduction 



Critical temperature is known only for very few spin models like the two-dimen- 
sional Ising model, which is exactly soluble in 2d, and Potts models where the 
critical temperature can be found by the duality relations [1-3]. In principle a 
position of the critical point can be obtained using one of many existing approximate 
techniques ranging from high and low temperature expansions to renormalization 
group methods. In this paper we employ a simple method based on the moments 
of the transfer matrix Q. Let T be the transfer matrix for a ci-dimensional spin 
system with linear size L and let f3 = l/Zc^T. One defines a characteristic function 
of rank n 



As was conjectured in [|| the maximum of p n should occur at the critical point 

Pmax ficrit- (2) 



The existing results suggest that the validity of the method is related to the dual 
properties of the system. It was checked directly for the Ising and Potts models with 
low q values on the square lattice that the criterion reproduces exact values of the 
critical temperature. It can be made exact also for the isotropic Ashkin- Teller model, 
for which the duality transformation involves three couplings. However, in this case 
the definition of the characteristic function should be modified [5,6]. On the other 
hand, for the two-layer Ising model, which is not self-dual, the criterion still gives 
surprisingly good numerical estimate of the critical temperature |7[]. Therefore it is 
very tempting to perform further investigations of the method to test its applicability 
in different situations. This is a part of our programme whose final goal is an exact 
criterion for not self-dual systems. 

In this paper we test the proposal (Eq. (|2])) for the two-dimensional Potts model. 
Contrary to the previous mainly numerical work || our discussion is based on the 
explicit analytical expression for the function p 2 for arbitrary q. We also check some 
consequences arising from the criterion in the limiting case n —>■ oo. These are 
interesting equalities relating internal energies of one and two dimensional systems. 
The calculation of p2 is given in Section 2 and the equalities mentioned above are 
derived in Section 3. Appendix contains some technical details. 



2 The p2 function 

The partition function of n coupled chains is related to the matrix T n , which prop- 
agates horizontally n spins 

Z n = TrT n = Tr(T n ) L , (3) 

where L is the linear extension of the system as in the formula ([]]). The construction 
of the matrix T n is shown, for n = 2, in Fig. 1 where periodic boundary conditions 
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are assumed in vertical direction. Since in the thermodynamic limit [L — > oo) only 
the largest eigenvalue, t njmax , of the matrix T n dominates, we obtain for p n the 
expression 

Pn((3) = f^, (4) 

which is very convenient for analytical calculations. 

For the Potts model the matrices T\ and T2 have a simple form and it is relatively 
easy to find analytical expressions for the largest eigenvalues. This is done in the 
Appendix where the formulas for the P2 function are given for square and triangular 
lattices. The functions P2(x), where x = e 13 , for the Potts model on the square lattice 
are shown in Fig. 2. For arbitrary q they have the same characteristic shape with the 
maximum, whose height and position depend on q. The position of the maximum 
is located at x — 1 + ^fq which agrees with the standard values obtained by duality 
of the six vertex model [1,3]. Indeed at this point the Potts model undergoes the 
second order transition for q < 4 and first order for q > 4. 

The triangular lattice is not self-dual but it is dual to the honeycomb lattice. 
The exact critical temperature still can be found from some properties of ice-type 
models and it is given as a root of the cubic equation [1,3] 

x 3 - 3x - q + 2 = 0. (5) 

The results for this case are summarized in Tab. 1, where we compare the critical 
values and the position of the maximum of p-i. The criterion still gives a good 
estimation of the critical temperature. In particular for q = 2 it is exact, which is 
related to the existence of the star-triangle relation which provides a certain kind of 
duality transformation. 



3 Critical equalities 

By taking the logarithm of (Eq. 1) and by differentiating it with respect to x one 
obtains the condition for the internal energies 

^l^max) Uni^X 'max) j 

where U\ (u n ) denotes the density of the internal energy for a spin chain (n coupled 
spin chains). In the limit (n — > 00) we conclude that at the critical point the energy 
density of the 2d system equals to the energy density of a single chain. This allows us 
to use standard expressions for the internal energy in Potts models to check validity 
of our criterion also in this limiting situation. 

Let us first consider the square lattice. The internal energy of a chain per one 
spin is: 

2x + q-l 

«i = — ; r- 7 

x + q-1 



(6) 
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Evaluating this at the critical point x cr u = 1 + y/q we obtain, in agreement with 
known results [1,3], namely that the energy of the 2d system is: 

Ucrit = 1 + —=■ (8) 

For q > 4 when the system undergoes the first order phase transition the internal 
energy at the critical point is discontinuous. To fix this ambiguity one usualy defines 
the energy as a mean of the high and low temperature phase, and so should be 
interpreted our results. 

In the last three rows of Tab. 1 we compare the predictions for numerical values 
of internal energies for the Potts system on the triangular lattice. E(T C ) and L 
are results for the internal energy and latent heat respectively taken from Ref. ||. 
E(T max ) is the energy of a single chain evaluated at the maximum point. The 
method still gives good numerical estimation of the internal energy. For the first 
order phase transition it fits nicely between energy of the high and low temperature 
phase. 

The critical equalities are simple mathematical consequences of the maximum 
rule but their physical origin is unclear. One can say that at the critical point when 
the correlation length £ diverges the interaction of a single chain with the rest of 
the 2d system can be simplified to the interaction of spins in a single chain with 
periodic boundary conditions. This would also mean that characteristic function 
should contain information about correlations in the system. Thus it should be 
possible to obtain other characteristics of the transition, for example to distinguish 
the order of the transition and to find the critical exponents. One can may be 
improve the method by combining it with the renormalization group techniques. It 
would also be tempting to look for signals of the order of transition containing in 
the function p n . For the Potts model shape of p2 exhibits rather slow crossover as a 
function of q and for the moment it is not clear weather one can see any particular 
change at q = 5. 
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4 Appendix 

In this Appendix we derive analytical formulas for the p2 function of the 2d Potts 
model. The energy of the q state Potts model is given by 

#00 = -£* w (9) 
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where the sum runs over nearest neighbours sites labeled by i, j and the spin variables 
s are assumed to take q distinct values. The transfer matrices can be specified by 
their elements between spin states. For the square lattice they are: 

(s|Ti|s') = exp/3(5 S)S / + 1), 
(sis 2 \T 2 \s' lS2 ') = exp(3(5 SuS2 + 8 sis[ + 8^ + 5^). (10) 

It can be readily verified that t\ max = x(x + q — 1). The largest eigenvalue of 
the matrix T 2 can be obtained noting the special structure of the Potts interaction. 
To make it more evident we introduce an index n = 0, . . . ,q 2 — 1, which enumerates 
successive rows (columns) of the matrix, n can be represented as a two digit number 
in base q. With each digit we associate a spin variable. Let us now construct the 
equation for the eigenvalues. The 0-th, (g + l)-th,. . . ,(q 2 — l)-th row of the resulting 
determinant we multiply by some constant, say A. At the end we add all the rows 
to the last one. It is easy to notice that we generically obtain two expressions in the 
last row: in the 0-th, (q + l)-th,. . . ,(q 2 — l)-th column w\ and in the other columns 
w 2 . This is because the Potts interaction depends on whether two spins are equal 
or different. Detailed form of W\ and w 2 can be obtained from some combinatorical 
considerations 

Wl = A(x 4 -X + (q- l)x 2 ) + (q- 2)(q - l)x + 2(q - l)x 2 , (11) 
w 2 = A{2x 2 + {q- 2)x) + x 2 - A + 2(q - 2)x + q 2 - 3q + 3. (12) 



If Wi = w 2 = one obtains a quadratic equation for a certain eigenvalue. We 
choose the larger root of this equation 
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Xi(x) = -(Q 2 (x) + ^Q 3 (x)), (13) 

where Q 2 and Q3 are polynomials given below. The other eigenvalues Aj(x) % = 
2,...,q 2 can be considered as continuous functions of x. Moreover Ai(l) = q 2 , 
Aj(l) = 0. By the Perron's theorem ||, X\(x) 7^ since the elements of T 2 

are positive. Thus Ai(x) must remain the largest eigenvalue in the whole interval 
(l,oo). 

Collecting our results together we finally obtain: 

P2{P) = 1 (14) 

where 

Qi(x) = 2x 2 (x + q-l) 2 , (15) 
Q 2 {x) = x 4 + qx 2 + (2q - A)x + q 2 - 3q + 3, (16) 
Q 3 ( x ) = x 8 + (2q - A)x & + (8 - Aq)x 5 + (-q 2 + 18q - U 
+ (l2q 2 - 32q + 16)x 3 + (2g 3 - 6q 2 - 2q + 12)x 2 
+ (4g 3 - 20g 2 + 36q - 24)x + (q 2 - 3q + 3) 2 . (17) 



)x 4 



The method described above applies also to the transfer matrices of the Potts 
model on the triangular lattice 



(s\T lS ) = ex P/ 3(2cW + l) 
( 8l 8 2 \T 2 \8 1 8 2 ') = expP(6 sl , S2 + 6 SiA + 6 S2S ,+6 sW2 + 5 Sisl2 + 6 S2s[ ). (18) 
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The final result has the same form as (|I4[) with 



Qi(x) = 2x 2 (x 2 + q- l) 2 , 

Q 2 {x) = x 6 + x 2 (q + 1) + x(4q -8) + q 2 - 5q + 6, 
Q 3 (x) = x 12 + (2q - 6)x 8 + (-8q + 16)x 7 + 



(19) 
(20) 



{-2q 2 + 26g - 28)x 6 + (17g 2 - 54g + 41)x 4 

+ (-8g 2 + 40g - 48)x 3 + (2g 3 + 12g 2 - 74g + 84)x 2 

+ (8g 3 - 56g 2 + 128g - 96)a; 

q A - I0q 3 + 37q 2 - 60g + 36. 



(21) 
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e /3cr 


v/3 


1.8794 


2 


2.1038 


2.1958 


2.2790 


2.3553 


4.8272 




V3 


1.8757 


1.9930 


2.0938 


2.1893 


2.2635 


2.3373 


4.7266 


E{T C ) 


0.8333 


0.7603 


0.7131 


0.6881 


0.6669 


0.6506 


0.6377 


* 


L{q) 











0.0310 


0.1172 


0.2042 


0.2795 


* 


E{T max ) 


0.8333 


0.7584 


0.7026 


0.6819 


0.6596 


0.6404 


0.6256 


0.4561 



Table 1: Comparison of the results taken from Ref. [3] and the values obtained 
from the maximum criterion for the 2d Potts model on the triangular lattice. 



Figure 1: Construction of the transfer matrix T2 for the two-dimensional Potts 
model on the square lattice. Periodic boundary conditions are assumed in vertical 
direction. 
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Figure 2: Functions pi for the Potts model on the square lattice. Three lines are 
shown for q = 2 (solid) , q = 3 (dashed) and q = 5 (dotted). 



